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1 Introduction 



The fragmentation of quarks and gluons into photons can be observed in e^e~ 
annihilation experiments and in the production of large-p^ photons in hadronic 
colhsions. This phenomenom is described by the distributions D'^{z, M^) and 
D2{z,M^) where z is the fractional momentum carried away by the photon and 
a time-like scale fixed by the hard process (M^ ~ in e+e~-annihilation 
and ~ p']_ for large-p_L photons). Unlike the fragmentation into hadrons 
which are complex bound-states, the photon has a known pointlike coupling to 
the quark. Therefore we expect these distributions to be fully calculable in per- 
turbative QCD. Witten was the first to show that this is indeed the case [1], at 
least for large enough to neglect non-perturbative effects, and he wrote the 
Leading Logarithm (LL) expressions for D2{z,M^) (a = q,g). In practice, it 
turns out that we need to know the fragmentation functions in kinematical do- 
mains where is not asymptotically large (M^ ~ ~ 25 GeV^ in fixed-target 
direct-photon experiments). These non-perturbative contributions and Beyond 
Leading Logarithm corrections (BLL) to Witten's LL results are sizeable. It is 
the purpose of this paper to present a careful analysis of these effects and a new 
parametrization of D'^{z, M^) and D^^z, M^). 

New experimental results justify the updating of an analysis published some 
years ago [2]. Since then the LEP collaborations studied in detail the quark 
fragmentation into isolated hard photons; the inclusive fragmentation functions 
D2{z,M^) should also be measurable [3]. On the other hand new precise data 
on direct photon production at large p± has been [4], or will [5] be presented 
soon, requiring more precise theoretical inputs. Finally it is now possible to 
better constrain the non-perturbative part of the fragmentation functions which 
is obtained by means of the Vector Meson Dominance Model (VDM). Indeed, 
new data on the inclusive p production at LEP p, |^ allow a better control of 
this contribution. 

The study of the fragmentation functions follows the theoretical approach 
developed in the analysis of the crossed reaction, namely the Deep Inelastic Scat- 
tering of a virtual photon on a real photon, which probes the parton contents of 
a real photon [8]. The good agreement between theory and data obtained in this 
channel lets us hope to derive sound predictions for the D1{z,M'^) distribution. 
Here we study only the inclusive fragmentation functions, without any isolation 
condition around the photon. The isolated case [9] raises theoretical problems 
concerning the Infra- Red stability of the prediction [10, 11] that we do not discuss 
in this paper. 
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2 Theoretical Background 

The fragmentation functions D2{z, M^) verify the inhomogeneous DGLAP (Dok- 
shitzer, Gribov, Lipatov, Altarisi, Parisi) equations [12, 13] (the convolution 
f ® g{z) is defined by du dv f{u) g{v)6{uv — z)) 



for the non-singlet sector {Cns,i = 2(e^— < >)), and 

Nf 

for the singlet sector (C, = 2Ar^ < >), with = + ^l) and L'^;^,^^ = 

i=i ' ' 

{D'^_ + D^.) — D'^/Nf. The inhomogeneous kernels have a perturbative expansion 



K^aiz.M 
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2^^^^°JW + ^(M^)i^Sn-) + ---) (3) 



as do the homogeneous kernels Pab- The kernels K.^a are given in [2], and the 
homogeneous ones can be obtained from [14]. Let us notice that the coupling of 
the gluon to the photon can only take place through a quark loop ; therefore the 
expansion (|^) of K^g starts at order 0{as). 

In the moment space (/(n) = dz z"-~^f{z)), (|l]) and (0) can easily be solved 
[15,16,17]. For instance for the non-singlet distribution we obtain 

D--(„. M^) ^ C ^^K,,(n) eJ-""' ^""^ (4) 



where we have dropped the suffixes ns and i. The suffix AN means anomalous, 
a qualifier given by Witten to the solutions of eq. in order to characterize 

their asymptotic behaviours. Indeed, with the definitions r = as{M'^) /as{MQ) 
and dn = 2P'^°)(n)//3o, the solution can be written, in the LL approximation, 

M') = -4%, — ^ (1 - r'-'-) (5) 



an expression which explicitly displays an asymptotic behavior proportional to 
In ^ (in (^ we kept only the lowest order term of the /5-function : 
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Expression (^) is not the full solution of the inhomogeneous equation (|l|) ; 
we can add to (§) a general solution of the homogeneous equation (eq. with 
Kq-y = 0), so that the full solution is 

D^n, M^) = M^) + M^). (6) 

The physical interpretation of expressions @) and (|^) is the following : D^{n, M^) 
is given by the sum of a perturbative component D^'^'^ and of a non-perturbative 
component D'^'^^ . D'''^^ is fully calculable in perturbative QCD, as long as 
is large enough, > Mq where Mq is the boundary between the perturbative 
and non-perturbative domain. For = Mg, the perturbative approach is no 
longer valid and is given by a non-perturbative fragmentation function D"''^^ , 
which verifies for > Mq the homogeneous DGLAP equations. 

The non-perturbative input D'^'^^{n, Mq) is not known. We modelize it fol- 
lowing VDM and we assume (for < Mq) that quarks and gluons first fragment 
into vector mesons which then turn into photons. Therefore we could write 

D^-^^(n,M2)=a C,D"{n,Ml) (7) 

where the fragmentation functions may be measured in e^e^-annihilation 
experiments. The coefficients are fixed by VDM. The value of Mg is not known, 
but should be of the order of an hadronic mass and we take Mq ^ ~ 0.5 GeV^. 
The same approach in the crossed channel 77* X leads to predictions in good 
agreement with data [8]. 

However the approach just described is too naive as it is based on a LL 
approximation. At BLL order, the decomposition is not factorization scheme 
invariant, and our VDM assumption (|^ must be refined. Let us study this 
problem which is related to BLL corrections to the LL expression (^. 



3 Non-Perturbative Input and Factorization 
Scheme 

We consider the one-photon inclusive cross-section in e"'"e~-annihilation. It is 
given by the convolution between the hard sub-process cross-sections Ca{z) and 
the parton fragmentation functions [18] 

Nf Nf 

° i=l i=l 

Nf 

+ 2j2etC^ (8) 

1=1 
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where (Tq = 47ra^/(5^. The hard cross-sections Cq{z) and Cg{z), which have 
expansions in powers of as{Q^) 

Ca = K, + ^iQ')ci'Hz) + --- (9) 

also appear in the one-hadron inclusive cross-sections and have been calculated 
in [19] (in (§) we consider the sum of the transverse and longitudinal cross- 
sections). C'^{z) is characteristic of reactions involving photons and describes the 
direct coupling of the photon to quarks in e"'"e~-annihilation. Its expression in 
the MS Factorization Scheme (FS) can be obtained from [19] : 

CJ = ^_n l^(ln(l-2) +21nz). (10) 

27r z 

Actually it is well known that the fragmentation functions and hard cross- 
sections are not univocally defined. For instance, a part of can be absorbed in 
the fragmentation functions, leading to a new photonic FS to which correspond 
new functions Da and C^. Each term on the rhs of (|^) is therefore FS dependent, 
but the sum is not, being a physical quantity. It is easy to verify that this implies 
that D'^''^^ in expression (^) is not FS invariant. Such ambiguities appear also 
in the definition of coefficients Cg and Cg and kernels Pat- The influence of this 
hadronic FS on the fragmentation functions was studied in ref [8]. In this article, 
we will focus only on the difficulties related to the photonic FS. 

In order to grasp this point more clearly, we calculate (^ including BLL 
corrections. Expanding in powers of a^, we obtain for expression (|^) 



2n A 




+ r-'^- D^^^^{n,Ml). (11) 

By combining this result with (^ and keeping the relevant terms proportional to 
(1 — r~°'") and r~'^", we easily obtain (for the non-singlet contribution and writing 
again Cns,i instead of C) 

1 da'^ , , o f ^ a 2 

In) ■■ 



(To dz 
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+ C„,, C-'(n) + r>,I;7(«, M^) r-*'} + 



a 2 



271 Po 



n] 



+ [Cns, C\n) + D::;7(n, Ml)] r~'- + ■■■} + ■■■ (12) 



Expression (|T2D explicitly shows [20] that the combinations C^{n)dn—^^K- 

and Cns,i C'^in) + D,^^^'^(n) are FS invariant ; if we change the scheme, we must 
obtain 

Cns, C^in) + &J/{n) = C^(n) + D^'^Yin) (13) 

and we clearly see that the "non-perturbative" component cannot correspond to 
a VDM contribution alone, which should be FS invariant. 

In ref [8], we discussed the structure of D"''^^ (actually in the DIS channel) 
in detail and showed that it consists of two parts. One part includes all the non 
perturbative effects and is scheme independent. The other part depends on the 
scheme and can be perturbatively calculated. It corresponds to the coUinear part 
of C^. In this paper we quote the result without proof, refering the interested 
reader to the original paper [8]. 

Dl^'/iz, Ml) = Di:f'\z, Ml) - C.., D^^%z) (14) 

where 

D^^{z) = — ( ^"^^^""^^^ (ln(l -z)+ \n{z)) + z] . (15) 
2n \ z J 



Expression ( jT^) is our initial condition at = Mq. A similar result can be 
obtained for the singlet sector with Cns,i replaced by Cs- For the gluon fragmen- 
tation, we have D]'^P{z,MI) = D]'^^^^ {z, Ml). 

The previous discussion is valid for light quarks. For massive quarks, we 
neglect the VDM component. For instance, we neglect the ?/;-dominance contri- 
bution to the fragmentation of a charm quark into photons. But we still have 
a "non-perturbative" input. Indeed Nason and Webber [21] calculated the frag- 
mentation of a heavy quark or anti-quark into a photon (actually a gluon in these 
calculations) with the result (dropping powers of rriQ / Q) 

where the direct term Cq, calculated in the massive FS, is given by : 

C^q{z) = C\z) - |- i±ii^(21n;^ + 1). (17) 
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C^{z) is the direct term in the MS scheme given by (|l3). By taking into account 
BLL corrections one obtains an expression similar to ([I^), but in which the kernels 
are calculated in the massive scheme and r is replaced by rq = as{M^)/as{mQ). 
In particular, in the terms 



-dn 



in) = -— ^K^^Mn) + CMn) + ■ ■ 



2e^ao dz ' 2n [5,^"^^"' d^ 

one recognizes the FS invariant combination 



a 2 <cj(n 
2vr (3q dn 



+ Cl{n). (19) 



Expression (|T8D can be transformed into the MS scheme we used in this paper : 

a 2 l-r""'" 



2n /3o ^ 4 



27r j3o\ ^"^^ ' ^ \ d, 




-dn 



a 2 6K{n)\ a 2 5K{n) 



with MS expressions 

-(1) _ 7^(1) 

and 



K« = K«+5if (21) 



From ( piD and ( P^ we see that we can recover the massive result (^Dj) by working 
in the MS scheme, but with a non-zero input at = tuq given by 

D-^(z, ml) = -£e^^^^^(2 In z + 1). (23) 

Let us end this section by comparing our present approach with previous BLL 
studies. 

The authors of ref [17] use a different approach but obtain similar results 
for the "non perturbative" input. Invoking the "perturbative stability", they 
choose to work with a factorization scheme (called DIS^) in which the direct 
term Cj)jg^{z) vanishes (more precisely the transverse direct term). Then they 
assume that the input at = Mq is given simply by the VDM contribution. 
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Using (|T3D, we can translate this input in the MS language. In the non-singlet 
case, we find 

p,^,NP _ f^-y,NP ^ r^l _ riiyDM ^ ^7 

where ~ now means DIS^ and with 

C:i = —^-^^^^—^(\n(l-z) + 2\nz)-2^—^ (25) 
2tt z z 

valid for any fiavour. We see that this expression is fairly similar to the one we 

obtained, namely eqs. (0) and (|23| ) ; it produces similiar effects when z goes to 

zero or one. 

The present approach differs from the BLL study of ref [2] in which the input 
^7, A/5 (^^^^2) is equal to zero. This leads to a different behavior of D'^{z, M^) at 
small values of z that we discuss in the next section. 



4 Numerical Studies of the Anomalous Compo- 
nent 

In this section, we perform a numerical study of the anomalous fragmentation 
function. Later we shall add the VDM contribution in order to obtain the com- 
plete fragmentation functions. We solve the DGLAP equation in which the kernels 
are massless and take into account the effect of the mass of heavy quarks by using 
thresholds at /i^ = and /x^ = ml with D2{fi'^ < ml) = and -D^(At^ < mf) = 0. 
Then, the input (|^) is introduced at m^ and mf,. We use the following values: 

^QCD = 230 MeV, mcharm = 1-5 GeV and mbottom = 4.5 GeV. 

In fig. 1 we display the anomalous fragmentation functions obtained with 
Mq = 0.5 GeV'^ and the input D'^'^^{z, M^), whereas in fig. 2 we show the 
results obtained with the boundary condition D'^'^^ {z, Mq) = i^D'^'^'^^ = 0). 
The effects of D"''^^^ are important at small values of z, especially in the gluon 
case. In both figures, the gluon fragmentation function is negative at small z. 
But the z-range in which the D^'^^ is positive is larger with the D'^'^^^ input. 

This small- 2; behavior of D'^{z, M^) is due to BLL corrections to the LL so- 
lution which does not show such a pattern. The BLL kernels have a singular 
behavior at small z : 

n'i'W - {cs(-4!i^).0(ili£)}. (26) 
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10^ I ^ . ^ . ^ . ^ . ^ 1 

0.0 0.2 0.4 0.6 0.8 1.0 
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Figure 1: Anomalous component with D^'^^ {z, Mq = 0.5 GeV'^) = D'^'^^{z) 
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Leading singularities in 1^^ removed 

Leading singularities in Wj", Pj" and Pj'^ removed 



Figure 3: Comparison of the anomalous gluon fragmentation functions with a null 
input at Ql = 0.5 GeV"^ and various singularities removed from kernels. 

The effect of the BLL inhomogeneous kernel K!y]] is particularly important, be- 
cause the Leading Order term vanishes {K^J =0). If we drop the most singular 
term (|26|) of K^^J (z), we obtain a gluon fragmentation function which becomes 
negative only at very small values of z {z < 3.10"^) (fig. 3) where the effect of 
the homogeneous kernels is important. When the singular behavior of Pgl^ and 
Pgl^ are also removed, the fragmentation function is positive. 

The z-domain in which the singular parts of the kernels are important has not 
been explored by experiment. At LEP, we have z > .7 and in large-p_L experiments 
< z >~ .5, far from the region where D^lzjM"^) is negative. Therefore it is not 
necessary to treat this small-z region more carefully by resumming to all orders 
the singular terms (0). 

When z — > 1, the kernels are also singular and the quark fragmentation func- 
tions are dominated by the BLL inhomogeneous contribution 

D^n M^) ~ -—e' ^^''^""^ — ^^'^/^ - (27) 
Ug{n,m ) 27r P(o)(n) 27r 2Cp In n 47r ^ n ^ ' 

showing that 

n- 1 

D-^{z,M^)r.-ln^^. (28) 
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In the cross section (H) this logarithmic term is cancelled by contribution coming 
from Cq^^ and ; as a result, the cross-section is regular when z — »• 1. 



5 Vector Dominance Model and 
Non-Perturbative Input 

In the Vector Dominance Model, the photon is described by a superposition of 
vector mesons (we neglect the J/ if) contribution) 

where ~ a. In e"'"e~-annihilation, the final quark (or antiquark) first fragments 
into a vector meson (or a [qq) state of spin 1) which is coupled to a photon through 
791). From (E9D we obtain 



where D^^^ is the fragmentation function of quark q into the (nn) bound state. 
We assume that the fragmentation of the quark q into the nn bound state is given 
by the fragmentation into a p-meson : 

jjim) = 2Df''' + Df'' 

Dfl = (31) 



DP ''" is the "valence" part for which the p°-meson contains the quark q and 
is the "sea" part for which the quark q does not enter the meson. The factor 2 
comes from the SU(3) wave function of the p°-meson. We can express the VDM 
fragmentation function (|30|) in terms of the quark and gluon fragmentation into 
p°-meson. 



We use data from ALEPH and HRS ^ {^/s = 29 GeV) in order to 



constrain Dp{z,M^) and DPiz.M"^). We found that data from MARK II 123 



TASSO 13 and DELPHI are not compatible with those from ALEPH and 
HRS. Because HRS has the greatest statistics, we have chosen the latter. Since 
data from JADE does not add constraints, they was not taken into account. 

We use the following parametrization of the fragmentation functions at — 
2 GeV^ for the gluon and the quarks up, down, strange and charm and at Qq = 
ml GeV^ for the quark bottom : 
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DP''{x) = D^/ix) = DP/{x) = Ns - xf'' 

DP''{x)= N^x^'^il-xf^ (32) 

D^'«(x) = NhX^^'il-xy^ 

DP/{x)= Ngx'^'il-xY' 



We reduced the number of free parameters in order to avoid too strong a 
correlation between them. We make the assumption that the behavior of the c 
and b quarks is related as follows : ai, = ac and [3h = /3c + 2. Furthermore, the 
exponents are fixed. First, we found that it is not possible to fit HRS and 
ALEPH data if we keep the HRS point at a; = 0.652. For this reason, we made the 
fits without this point. Then, in a first fit, we fixed jSg (set I). When comparing 
the ratios D'^^ / D'^ (we used results from ref [^]) and D^/Dp^, we noticed that 
the former is bigger by a factor of 3 to 10 (depending on the value of x and Q"^) 
than the latter. Because the difference between non-perturbative mechanisms of 
fragmentation into or 7r° should be reduced in these ratios, they should be of 
the same order. Therefore our gluon fragmentation function which is not well 
constrained by our e"'"e~ data is probably too small. Thus we performed a second 
fit for which we fixed Ng in order to obtain a ratio D^/Dp of the order of the 
same ratio for pion (set II). The values of the parameters are shown in Table |l]. 
We can see that the increase of Ng implies a decrease of the normalization for 
the heavy quarks. We plot on figures ^ and |^ the comparison between fitted data 
and computed cross sections. The Xdof equal to 1.33 for set I and 1.22 for set 
II. 



set I 


set II 




N 


a 


P 




N 


a 




valence (u,d) 


0.785 


-0.5 


1.499 


valence (u,d) 


1.140 


-0.2 


1.693 


sea (u,d,s) 


0.111 


-1 


2.912 


sea (u,d,s) 


0.100 


-0.3 


3.437 


c 


0.567 


-1 


5.502 


c 


0.132 


-1 


4.820 


b 


1.020 


-1 


7.502 


b 


0.103 


-1 


6.820 


g 


0.108 


-1 


3 


g 


2.550 


-0.3 


3 



Table 1: Fitted parameters for the fragmentation functions to p. The exponents 
tta are fixed. 
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0.0 0.2 0.4 0.6 0.8 1.0 

Figure 4: Comparison between ALEPH data and predictions corresponding to set 
I and II. Black dots correspond to points used in the fits. 




0.0 0.2 0.4 0.6 0.8 1.0 



Figure 5: Comparison between HRS data and predictions corresponding to set I 
and II. Black dots correspond to points used in the fits. 
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6 Full Fragmentation Functions and Compari- 
son with Experiment 

We obtain the complete fragmentation functions by adding the VDM contribu- 
tions to the anomalous contributions. They are given in fig. ^ ^ (for = 
100 GeV^) and fig. |, | (for = 10^ GeV^). 

The results of fig. ^ and |^ correspond to set I of parton into p-meson frag- 
mentation functions discussed in the preceeding section, whereas those of fig. |^ 
and 1^ correspond to set II. We notice a sizeable difference only for the gluon frag- 
mentation functions; in this case the VDM contributions are very different. On 
the other hand the VDM contributions to the quark fragmentation functions are 
small and the curves of fig. ^ to ^ are very similar to the corresponding curves of 
fig. These distributions can be compared with those we obtained in ref [2] (fig. 
10 to 13). The latter are larger at small z and not too large Q^, the difference 
being essentially due to a different VDM input. In ref [2], we assume that the 
fragmentation in p-meson is similar to the fragmentation in vr" and we use the 
distributions D'f of ref [7] as VDM input. In this case, the input we obtained 
after a fit to data is much smaller. 

In fig. |iy, we compare our results with the LL parametrization of Duke and 
Owens [12]. However one must keep in mind that BLL distribution functions are 
factorization scheme dependent, and that our distributions are calculated in the 
MS scheme. A better comparison is provided by the cross section da''/dz, an 
invariant observable which can be compared to experiment. 

At present, there is no data with which to compare. ALEPH data could 
seem to be a basis for comparison, but it is produced by an analysis in jets. 
Following [3], this Collaboration defines the fragmentation into a photon within 
a jet. The fragmentation function Dj^^ defined in that way does not correspond 
to the functions calculated in this paper which are fully inclusive ; we do not put 
any limitation on the phase-space of the hadrons which accompany the photon. 

In order to better understand the difference between ALEPH results and our 
predictions, let us consider the decay of a Z-boson (of momentum Q) into a 
photon (pi), a quark (^2) and an anti-quark (ps). We define Zi = 2pi.Q/Q^ and 
yij = 2pi.pj/Q^, where Zi is the inclusive photon fragmentation variable. We have 
1 — Zi = (i,j,k different) and Ylii<jyij — 1- The variable used by ALEPH to 
describe the photon in the jet (here the jet is made of the photon and the quark) 
is = zi/{zi + Z2) = Zi/{1 + 1/12). However an integration is performed on yi2 
within the jet, so that the effective value zl-^^ at which one should compare our 
results is larger than z^. But if we assume that the largest contribution to the 
integral comes from the coUinear region yi2 ~ 0, we obtain zi ^ z^. Hereafter we 
use this assumption. 

One must also notice that the ALEPH Collaboration uses the Durham algo- 
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^-^^^ gluon 









0.0 0.2 0.4 0.6 0.8 1.0 
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Figure 10: The fragmentation function for gluon and quark up of our sets of 
fragmentation functions compared with those of Duke- Owens at = 100 GeV^. 
For the quark, set I and set II cannot be distinguished. 



rithm [28] to define a jet. According to this algorithm, y^°'^ = (1 — z^)/{l + z^); 
therefore the scale in the fragmentation function is no longer Q^, but (1 — 
z^)/(l + z^)Q'^ / (the extra 1/ z^ comes from the fact that in the inclusive case 
^max _ ziQ'^ and this factor zi ~ z^ is already included in our calculation). 

Finally, one must keep in mind that higher order QCD corrections to the qq^ 
process can generate logarithms of the jets parameter (e.g. Ucut) coming from the 
limited phase space integration, which are not present in the fully inclusive case. 
However for z^ large enough, Ucut does no longer constrain the phase space and 
the comparison between our predictions and ALEPH results is not spoiled by this 
effect. 



This comparison is shown in fig. [TTjfor 2-jets events and ycut = 0.1. We display 
our predictions for two scales in order to exhibit their sensitivity to the latter. The 
agreement is quite satisfactory. It is interesting to notice that, in this 2;^-region, 
one essentially tests the anomalous component of the fragmentation functions; 
once Qq is chosen, these parts are a pure prediction of the perturbative QCD. 
Qg, which is of the order of m^, caracterizes the border between the perturbative 
and non perturbative regions (cf eq.(6)). 
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flp=I/z^(l-z/(l+z.,)Q^ 



0.60 0.65 0.70 0.75 0.80 0.85 0.90 0.95 1.00 



Figure 11: Comparison between ALEPH data and our predictions for the direct 
production of photon. 



7 Conclusions 

We have studied the parton to photon fragmentation functions beyond the leading 
order. We recalled that the traditional decomposition of this functions in non- 
perturbative and anomalous parts depends on the photonic factorization scheme. 
Performing a careful analysis of this dependence, we propose a new definition 
of the perturbative and non perturbative components. In this approach, all the 
scheme dependence is put in the perturbative part. By using a VDM approach, 
we constrain the non perturbative component of the fragmentation function that 
we deduced from parton to rho fragmentation functions. The latter was obtained 
from a fit to LEP and PEP data. Finally, we propose two new sets of parton to 
photon fragmentation functions We used them to give new predictions for the 
production of direct 7 at LEP that agree well with experimental data. However, 
these data obtained by an analysis in jets are not fully inclusive and they do 
not exactly coincide with our inclusive fragmentation functions. Therefore, fully 
inclusive data for direct photon production would be very interesting, as they 
would allow to test a beautiful prediction of perturbative QCD. 



^FORTRAN subroutines which compute these fragmentation functions are available on re- 
quest by e-mail to fontanna@qcd.th.u-psud.fr 



18 



References 



[1] E.Witten, Nucl. Phys. 210 (1977) 189. Actually Witten obtained his results 
for the crossed reaction, the Deep Inelastic 77* Scattering, by using the 
operator-product-expansion method. 

[2] P. Aurenche, P. Chiappetta, M.Fontannaz, J. Ph. Guillet and E. Pilon, Nucl. 
Phys. B399 (1993) 34. 

[3] D. Buskulic et al., ALEPH Collaboration, Z. Phys. C69 (1996) 365. 

[4] E706 Collaboration, T. Fcrbcl in "Proceedings of the Rencontres de 
Moriond" , March 1996, edited by J.Tran Thanh Van. 

[5] UA6 Collaboration, private communication from M.Werlen. 

[6] P. Abreu et al., DELPHI Collaboration, Z. Phys. C65 (1995) 587. 

[7] D. Buscuhc et al, ALEPH Collaboration, Z. Phys. C69 (1996) 379. 

[8] P. Aurenche, M.Fontannaz, J.-Ph. Guillet, Z.Phys. C64 (1994) 621. 

[9] G. Kramer and B. Lampe, Phys. Lett. B269 (1991) 401 

E.W.H. Glover and W.J. Stirling, Phys. Lett. B295 (1992) 128 
Z.Kunst and Z. Troczanyi, Nucl. Phys. B394 (1993) 139 

[10] E.L. Berger, X. Guo and J. Qiu, Phys. Rev. Lett. 76 (1996) 2234 

[11] P. Aurenche, M. Fontannaz, J. P. Guillet, A. Kotikov, E. Pilon, Phys. Rev. 
D55 (1997) R 1124 

[12] C. H. Llewellyn Smith, Phys. Lett. 79B (1977) 189. 

K. KoUer, T. F. Walsh and P. M. Zerwas, Z. Phys. C2 (1979) 197. 

A. Nicolaidis, Nucl. Phys. B163 (1980) 156. 

D. W. Duke and J. F. Owens, Phys. Rev. D26 (1982) 1600. 

[13] We follow the notations of [2] . 

[14] G. Curci, W. Furmanski and R. Petronzio, Phys. Lett. B175 (1980) 27. 
W. Furmanski and R. Petronzio, Phys. Lett. B97 (1980) 437. 

[15] W. Furmanski and R. Petronzio, Z. Phys. C (198 ). 

[16] In this paper, we use a FORTRAN program written by P. Nason and modi- 
fied to take into account the inhomogeneous kernels. 



19 



[17] M. Gliick, E. Reya and A. Vogt, Phys. Rev. D48 (1993) 116. 
The moments of the kernels can be found in this reference. 

[18] W. A. Bardeen and A. J. Buras, Phys. Rev. D20 (1979) 166 derived expres- 
sion (8) for the DIS case. 

[19] G. AltareUi, R. K. EUis, G. MartineUi, S.Y. Pi, Nucl. Phys. 160 (1979) 301. 

[20] An all order discussion of the FS dependence can be found in appendix B of 
[8]. 

[21] P. Nason, B. R. Webber, Nucl. Phys. B421 (1994) 473. 

[22] S. Abachi et al., HRS Collaboration, Phys. Rev. D40 (1989) 706. 

[23] H. Schellman, MARKII Collaboration, SLAC-PUB-3448 (1984). 

[24] W. Brandehk et al, TASSO Collaboration, Phys. Lett. 117B (1982) 135. 

[25] W. Bartel et al., JADE Collaboration, Phys. Lett. 145B (1984) 441. 

[26] J. Binnewies, B.A. Kniehl, G. Kramer, Phys. Rev. D52(1995) 4947 

[27] P. Chiappetta, M. Greco, J.P. Guillet, S. RoUi, M. Werlen, Nucl. Phys. B412 
(1994) 3 

[28] Y.Dokshitzer, Z.Phys. C17 (1991) 1441 



20 



